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In  studying  of  flows  in  the  region  of  bran-sonic  speeds,  one  has  to  study  prob- 
lems described  by  non-linear  differential  equations.  In  connection  with  this,  basic 
results  in  solving  such  problems  were  obtained  by  a time  curve  method  and  a method 
of  series  approximations.  The  difficulties  which  appear  using  these  methods  are  well 
known.  A nWnh  f~r  general  methods  of  solution  led  to  development  of  series  of  new  appro- 
ximation methods  of  study  of  own  country  problems  for  transonic  flows  (1-3). 

The  asymptotic  method  studying  problems  of  around-sonic  non-accented  gas  dyn- 
amics, based  on  a problem  of  full  approximation  (6)  was  developed  in  the  works  (^,5)» 
Applying  his  method  to  a series  of  non-linear  problems  showed  its  effectiveness,  co- 
munity  and  reliability  of  results(4,5,7) - 

Further  development  of  a complete  approximation  method  was  given  in  the  said  work, 
and  on  its  basis  was  shown  a possibility  to  investigate  non-linear  problems  of  superso- 
nic  range  of  supersonic  flows  at  M^  > 1,  such  that  M^-  1 = 0 (fc),  fc<£-l  - 

small  parameter  of  destruction.  The  method  has  a sufficient  community  of  aim,  which 
makes  it  possible  to  study  space,  axis  symmetrical,  conical  and  other  problems  of 
sonic  flows. 

In  applying  the  asymptotic  methods  to  mechanical  problems,  important  th§  questien 
of  real  possibility  to  construct  algorithms  of  higher  approximations.  The  basic  comple- 
xities of  calculation  of  non-linear  effects  are  connected  with  a great  labour-consuming 
realization  of  higher  algorithms  approximations,  so  that  a problem  to  construct  appro- 
ximation above  the  second  is  often  hopeless.  On  the  other  hand,  corresponding  final  in- 
formation of  exactness  of  the  mathematical  model  in  a series  of  problems  determines 
the  contentment  of  second  approximation.  So,  the  er.trophy  change  in  a transonic 
flow  is  in  the  order  of  0 ( e*’)  and  the  assumption  ol  flow  isotrophy  determines  the 
contentment  of  second  approximation.  Simultaneously,  the  construction  of  a third  appro- 
ximation requires  entrophy  calculation,  which  significantly  complicates  problem  solving. 

In  connection  with  this  , there  is  a current  problem  to  construct  effective  algo- 
rithm and  methods  of  construction  of  second  approximation.  With  regard  to  this,  the 
problem  of  full  approximation  was  set  as  a problem  to  select  a physical  space  roliec- 
tio . in  some  new  space,  in  which  the  solution  of  the  first  approsimation  would  give, 
at  a transition  into  a physical  space, solution  of  two  first  approximations.  The  reflexion 
is  realized  by  means  of  a special  deformation  of  coordinates,  which  transforms  a non- 
linear differential  equation  in  a physical  space  into  a linear  in  auxiliary  space,  a 
space  of  approximation.  Therefore,  the  method  of  full  approximation  is  by  ideology 
close  to  the  method  od  deformated  coordinates  (8).  The  central  question  in  a problem 
of  full  approximation  is  a proof  of  its  solution  for  a type  of  equation  which  is 


I 


being  studied.  In  this  work, problem  solving  of  full  approximation  of  non-linear  dif- 
ferential equation  of  supersonic  gas  dynamics  in  a physical  space  is  investigated  by 
means  of  linear  equation  of  space  approximation  which  is  proved  by  the  following  theo- 


Theorem  on  asymptotic  relation.  Let  a large  number  of  elements  E (^2)  G C ( Q ) 
satisfy  the  equation  (9) : 

+ V - (K  - 1 ) <P„  = Ml  [ ^ (2tpx  -f  4-  rl  + tfx 

x <»«  + V + ^ + (2V*  + Vp9.cz  + vlvn  + 1 pl<pu  + (]) 

2W„  + 2(‘+  Vz)V,9„  + 2(1+  <px)  ] ; 

?€ft;  V£E( Q); 

and  number  of  elements  E (5c*)  £ (Q()  satisfy  the  equation 


••  • 


ftc  + V„-9tt  + 0(e2+a)  =0;  *,€!!,;  9€£,<Q1). 


If:  l)  a number  of  elements  E ((S^.)  and  E^  (^?()  are  of  eaual  power; 

% NL  *+  t ^o  ve);  X = 0 (e^2);  v = 0 (e0/l);  0 < a < 1 ; 

3>  9 = «p;  H-vll — 1;  S€U  0 < n < 1; 

4)  operator  P,  which  realizes  the  reflection,  gives  the  relation  of  elements 

<ftx,  y,z)  i <7,  (I.  *1.  t) « in  the  form 


p + 0;  r#-‘;  »-'}) 


fr-Xg;  P = (A&-1),/2; 

This  has  an  asymptotic  connectio  between  the  solutions 

9 (x,  y,  z)  = 9 (I,  rj.  S)  + 0 (eJ+°), 

0 

where:  0(e2+a)  has  characteristic  l>m  — = — -N, 

e-*0  e 

N is  certain  constant 
Q i Qi  are  regions  of  the  potential  flow 
q qnd  q^  are  regional  points 
e are  dimentionless  parameters 

a,  n are  exponents,  certain  from  condition  of  limits. 

8^(1.  „ * . . 

Result  No.l.  For  small  distinctions  * and  the  transition  (.3)  transfers 

into  transition  of  Prandtl  - Glauert,  and  the  equation  (l)  transfers  into  a linear 
differential  equation  of  small  destructions  theory. 

Result  No. 2.  In  selecting  of  transition  (3)  in  the  form  of  Prandtl  - Glauert, 
the  equation  (l)  transfers  into  a non-linear  differential  equation  of  small  destruction 
theory  for  transonic  flow  (2). 

[A4*  — 1 + (V  + »)  <P*I  <Px*  — — V*.  + 0 *=  p* 


o 
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The  proof  of  the  theorem  is  analogous  to  the  proof  of  the  theorem  for  non-accetuated 
flows  (4,  5) and  is  not  shown  here. 

As  an  example,  let  us  investigate  a flow  problem  of  a bodily  profile, at  an  angle 
of  atack  and  Mach  value  a little  larger  than  unity, by  an  isotropic  transonic  gas  flow. 
The  angle  of  Mach  number  and  the  profile  edge  are  such  that  on  the  front  and 

rear  edges  are  formed  joined  jumps  of  condensation  with  a supersonic  flow  at  both  sides 
of  jumps  and  the  deformation  function  B^-,  (+  1,  0,  0)  = 0.  In  this  case,  the  conditions 
in  the  theory  of  the  second  approximation  are  not  used  on  jumps  in  a visible  form.  If 
the  form  equation  is y*  =6  T (x )\~F  = 0(!),  the  problem  in  a physical  i..yi.  a’-ea  will  be 

l<-  1 +A*L(Y+  1)9.1^-^  +0(e2+“)  = 0;  q€Q; 

=(1  + 9**)^,.;  x£\— 1.-9-+  1);  q€f  ; (5) 

9 = 9 =0;  |x|>  1. 

Using  the  evaluations  2,  3 from  the  theorem  of  a limiting  condition,  we  shall 

2 

obtain  an  expression  for  exponent  n = -g  at  which  a = if  shall  have  known 

value  of  n-2/3  (10).  Since  the  flows  m an  upper  (+)  and  lower  (-)  areas  are  indepen- 

dent, let  us  investigate  the  upper  area.  According to  the  characteristic  equation  theo- 
rem in  the  area  of  approximation  in  physical  area  is 

^ dy  I (,  \ 

Tx  3(ML-1)  v 

Integrating  this  expression,  we  shall  obtain  equation  of  curves  if  = const 


* 3(AC-A):.V 

We  shall  write  the  solution  of  the  problem  equation  (5)  in  the  form 


H±JLi  \ 


Using  limiting  condition,  we  shall  obtain  and  expression  for  <^that  means  also  for 
the  pressure  coefficient: 


+ 1 


(An*—  i) 


Thefj  results  entirely  'egrae  with  the  known  results  of  the  second  approximation 
theory  of  Buseman  and  Van-Ih/ which  were  obtained  by  method  of  consecutive  approxi- 
mations (9)  and  method  of  deformed  coordinates  (8).  However,  our  solution  has  a consi- 
derable overweight;  it  is  evenly  useful  in  all  of  the  flow  area,  providing  that  it  does 
not  contain  parts  which  lead  to  a cumulative  effect  (8)f  and  at  the  same  time  gives  an 
equation  of  corrected  Mach  lines. 

It  should  be  pointed  out  that  all  qualitative  results  , obtained  in  a no:i-'.inear 
1j.ow  of  profiles,  come  from  the  obtained  formulas  (11).  The  obtaining  of  the  non-li- 
near result  already  in  the  solution  of  the  first  approximation  makes  the  method  of  full 
approximation  approach  the  method  of  accelerated  convergence  (12).  The  proof  of  the 
theorems  (4,  5)  together  with  the  one  given  in  this  paper  is  constructive;  it  gives 
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also  a selection  of  presentation.  To  the  discussed  solution  , one  may  apply  also 
the  method  of  the  ( 4,  5»  7 ) works.  The  procedure  of  solving  in  this  case  is  some- 
what longer,  however,  it  leads  to  the  same  results. 

The  overweight  of  the  method  that  is  developed  is  apparent  if  one  compares  the 
method  of  full  approximation  with  the  method  of  accelerated  convergence.  In  the  method 
of  accelerated  convergence,  beginning  with  the  second  step  one  has  to  find  integrals 
of  Pausson  equations  or  nonhomogenous  wave  equations,  which  is  connected  with  great 
calculating  difficulties  (9,  13)  and  , in  addition,  it  leads  to  non-homogenous  solu- 
tions, connected  with  * cumulative  effect.  Formalization  of  the  full ‘‘^fp ro x imat ion  me  - 
thod  by  itself  already  determines  the  problem  if  liquidation  of  these  integrals.  This 
significant  overweight  of  the  method,  at  last,  determines  its  great  possibilities. 
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